Oscillatory ac- and photoconductivity of a 2D electron gas: Quasiclassical transport 

beyond the Boltzmann equation 
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We have analyzed the quasiclassical mechanism of magnetooscillations in the ac- and photocon- 
ductivity, related to non-Markovian dynamics of disorder-induced electron scattering. While the 
magnetooscillations in the photoconductivity are found to be weak, the effect manifests itself much 
more strongly in the ac conductivity, where it may easily dominate over the oscillations due to the 
Landau quantization. We argue that the damping of the oscillatory photoconductivity provides a re- 
liable method of measuring the homogeneous broadening of Landau levels (single-particle scattering 
rate) in high-mobility structures. 
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I. INTRODUCTION 

An intriguing development in the study of a high- 
mobility two-dimensional electron gas (2DEG) was the 
recent observationi of magnetooscillations of the photo- 
conductivity of a sample subjected to microwave radia- 
tion, as a function of the ratio uj/lUc- Here uj and lUc are 
the radiation frequency and the cyclotron frequency, re- 
spectively. Subsequent experimentS)2*^i^ working with 
very-high- mobility samples, yielded yet another dramatic 
discovery: for sufficiently high radiation power, the 
minima of the oscillations evolve into "zero-resistance 
states", i.e., the dissipative resistance of a sample be- 
comes vanishingly small. 

The nature of the oscillations in the photoconductivity 
(Tph has raised a lot of interest. An important step was 
made in Ref. 0, where a direct connection between the 
emergence of the zero-resistance states and the oscilla- 
tions was emphasized. Specifically, it was recognized that 
whatever the nature of the oscillations, when they become 
so large that the linear dc response theory yields a nega- 
tive fjph, an instability is developed leading to the forma- 
tion of domains of counter-flowing currents and thus to 
the zero measured resistance. Following this approach, a 
key issue which needs to be settled for understanding the 
experiments is the microscopic mechanism of the oscilla- 
tory photoconductivity (OPC). 

Similarly to the conventional Shubnikov-de Haas oscil- 
lations, the growing body of theoretical work is focused 
on the oscillations of the density of states (DOS) induced 
by the Landau quantization as an essential element of 
the construction. The mechanism of the OPC identi- 
fied in Ref. and analyzed in detail in Ref. |3 hinges on 
the oscillations of the DOS and is related to a radiation- 
induced change of the electron distribution function in 
energy space, /(e), such that /(e) oscillates with vary- 
ing both s/uJc and uj/uJc- A hallmark of this contribution 
to (Tph is that it yields an amplitude of the OPC which 
is proportional to the inelastic relaxation time Too due to 
electron-electron collisions (more effective at low temper- 



atures than electron-phonon scattering). Another mech- 
anism of the OPC, based on the effect of radiation on 
impurity scattering in the presence of the Landau quanti- 
zation, was suggested in Ref. (an earlier, closely related 
variant of this approach was formulated in Ref. llOj) . A 
systematical theory of this contribution to dph was con- 
structed in Ref. El Comparing the results of Refs. 00 
and Ref. Ull one sees that the mechanism^i^ dominates, 
i.e., leads to much stronger oscillations, if Tee ^ Tq, where 
Tq is the single-particle relaxation time due to impurity 
scattering. For typical experimental parameters, a char- 
acteristic ratio Too/ Tq ~ lO^jii Overall the results of ReL|^ 
are in good agreement with the experimental data as 
regards the behavior of cTph in the range of parameters 
where the OPC is not too strongly damped, i.e., where 
the experimental efforts have been focused so far. In 
particular, Ref. 8 explains the emergence of strong oscil- 
lations and, in combination with Ref. 0, the formation of 
zero-resistance regions. 

While the agreement between theory and experiment 
is very encouraging, the situation is not so clear in the 
experimental limit of weak (strongly damped) oscilla- 
tions. Central to the identification of the microscopic 
mechanism of the oscillations is, on top of their period 
and phase, the behavior of their envelope with decreas- 
ing magnetic field B. For any mechanism based on the 
DOS oscillations, the relation between the DOS and OPC 
damping factors is critically important. The OPOS'S'ii is 
damped at low i? by a factor exp(— Tr/ijJcTph), where the 
ratio Tph/Tq = 1/2 is a distinctive feature not sensitive 
to microscopic details of either disorder or weak inelas- 
tic interactionsiii^ However, as emphasized in Ref. 0, the 
experimentally reported values of Tq and Tph do not sat- 
isfy this relation, with Tph noticeably larger than Tq/2, 
roughly by a factor of 10 in Ref. and by a factor of 3 in 
ReL 2. Taken at face value, the difference would mean 
that the amplitude of the OPC observed at small B is 
orders of magnitude higher than given by the mechanism 
based on the Landau quantization, which might be con- 
sidered as a hint about a different origin of the OPC at 
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small B. Alternatively, the experiments on the damping 
of Shubnikov-de Haas oscillations might overestimate the 
single-particle scattering rate r~^, e.g., because of inho- 
mogeneous (due to macroscopic inhomogeneities) broad- 
ening of Landau levels. To resolve this dilemma, it is 
desirable to examine a range of mechanisms of the OPC 
in the absence of the DOS oscillations. 

In this paper, we analyze a mechanism of the OPC gov- 
erned by quasiclassical memory effects. These are related 
to non-Markovian correlations in electron dynamics^ 
We assume that jr/uJcTq » 1 and completely neglect 
weak oscillations of the DOS. The OPC induced by the 
memory effects is not specific to any particular type of 
disorder; however, below we concentrate on the follow- 
ing two-component model, "'^^ where the memory effects 
are particularly prominent. We assume that there is a 
smooth random potential of remote donors that are sep- 
arated by a large spacer d':$> kp^, where kp is the Fermi 
wavevector, from the 2DEG plane and, in addition, there 
are rare short-range scatterers, e.g., residual impurities 
located at or near the interface. We consider the case 
Ts ^ tl, where rg and tl are the zero-B momentum 
relaxation times due to the short-range scatterers and 
the long-range disorder, respectively. From the exper- 
imental point of view, this choice is motivated by re- 
ports (see, e.g., Ref. 15') that the zero-B mobility in very- 
high-mobility structures is frequently limited by residual 
impurities and tl/ts can be as large as 10. Although 
Ts ^ Tl in our model, we assume that Tq is determined 
by the smooth disorder, i.e., Tq ~ Ti/(2fc_F(i)^ <C ts- 

The paper is organized as follows. First, in Sec. m 
we outline the approach to the photoconductivity based 
on the Boltzmann equation. In Sec. IIIII we discuss the 
mechanism of the photoconductivity related to electron- 
electron interactions. In Sec. IIVI we turn to the magne- 
tooscillations induced by the memory effects. Our cen- 
tral results are presented in Sees. IVIVTI Section Ivl deals 
with the oscillations in the ac conductivity. Finally, in 
Sec. IVII we compare two mechanisms of the oscillatory 
photoconductivity, quasiclassical and quantum, related 
to the memory effects and the Landau quantization, re- 
spectively. 



II. PHOTOCONDUCTIVITY: ESSENTIALS 

A necessary input to the calculation of the quasiclas- 
sical OPC is the memory effects, discarded in the Boltz- 
mann equation. However, to set up a systematic formal- 
ism, it is instructive to begin with a derivation of dph 
within the conventional kinetic theory. The Boltzmann 
equation for the distribution function g(j), (f), t) of elec- 
trons in momentum space reads: 



Lg{p,cj),t) = -¥ dpg{p,(j),t) 



(1) 



where L = dt+tOcdj,- Id - /in, F 
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FIG. 1: Graphic representation of various contributions to 
the photoconductivity [Eqs. in the {v,n) space at the 

lowest order Uph ~ 0{£^^£^). 



of the momentum p with respect to the direction of £dc, 
/oi and /;„ are the elastic and inelastic collision integrals, 
respectively. 

We expand the distribution function at energy e in a 
series: g(p,4>,t) — ^^^^ (7^„(e) exp(w0 -f mwi). Elastic 
collisions lead to relaxation of angular harmonics with 
0; in particular, Icigui = —T^^gin, where t is the 
momentum relaxation time. Inelastic electron-electron 
collisions tend to equilibrate electrons among themselves 
but are not capable of establishing a steady-state dc pho- 
toconductivity. For the quasiclassical OPC (in contrast 
to that based on the DOS oscillations, cf. Ref. Q), the 
inelastic transitions due to electron-electron interaction 
do not play any essential role and will be neglected. To 
dissipate energy absorbed from the ac field, we intro- 
duce coupling to a thermal bath, e.g., to an equilibrium 
phonon system, characterized by a relaxation time Tin. 
Under the assumption that both Tin and the momentum 
relaxation time due to the coupling to the bath are much 
longer than r, the main role of the inelastic scattering is 
to yield a slow relaxation of the isotropic [v — 0) part 
of g to the equilibrium Fermi distribution fp at a bath 
temperature T. 

Expanding the nonequilibrium distribution function in 
powers of the driving force, we have 



g 



-L-^Fd^rfF 



(2) 



£dc is the dc electric field, £^ is the ac field, (j) is the angle 



A useful way of visualizing this solution of Eq. ^ at 
given order in S^c and £^ is by counting all possible 
couplings of harmonics g,yn represented as points on the 
(j/, n) plane (Fig. The dc field fdc couples nearest- 
neighbor harmonics along the v axis, gi/n gv±\.n- 

We 

are interested here in the linear (with respect to fdc) 
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photoconductivity Cph, so that only one such Hnk is al- 
lowed. The ac field £^ couples harmonics along diago- 
nals, gyn ^ gu±i,n±i and 5,.±i,„^i. The propagator L^^ 
is a diagonal matrix in (z^, n) space. The static longitu- 
dinal current j — cTph^dc, 



depRegio , 



(3) 



is expressed through gio, i.e., is given by a sum of all 
paths starting at (0,0) and ending at (1,0). Already at 
order a^h ~ C('^dc^w) many as 9 different graphs arise, 
shown in Fig. ^ to which one should add their counter- 
parts mirrored in the horizontal axis, which corresponds 
to the change lo — s- —lu. 

Our strategy for finding gio is to select graphs involving 
couplings whose strength diverges at rin/r — *■ oo. These 
are graphs returning to the point (0, 0) [graphs (a), (b) in 
Fig. which are proportional to (L^^)oo. This means 
that to order £^ the path (0, 0) (1, 1) ~^ (0, 0) ^ (1, 0) 
(and its counterparts in other quadrants) gives the main 
contribution to (Tph for 



Tin > r 



(4) 



The perturbative expansion in powers of £^ proceeds 
by iterating the loop (0,0) -> (1,1) -> (0,0). In this 
way we arrive at a simple relation (iiwc + '''^^)9±i.q = 
efdc dpf /2, where / = goo satisfies the closed equation 



2m 



(5) 



The function — + describes the absorption 
rate at energy e: 



approximation which neglects the e dependence of r in 
the integrand of Eq. ^ [normally, T(e) changes on a scale 
of ep] yields an identically zero photoresponse. Indeed, 
in that case Cph is equal to the static Drude conductiv- 
ity independently of the detailed shape of /(e), since 
J deedef = —2Trh?ne/m due to particle number conser- 
vation (rie is the electron concentration). It follows that 
the dependence of r on e should be taken into account. 
It is worth mentioning that, contrary to a naive expec- 
tation, this does not lead to any additional smallness of 
CTph since T(e) enters the result through expressions of 
the type e9eT|e=ep r. 

We do not discuss specific microscopic models of the 
inelastic coupling of electrons to a thermal bath, our pur- 
pose here is to use the simplest possible representation 
of /in- In a conserving relaxation-time approximation 



\\f-fF) , 



(10) 



where the e independent t-^^ is in general a functional of 
/(e), we get from Eqs. (0,13 for Te < ep: 



D ^ '^uj'm j^ii 



(11) 



Here K'^ 



dlKo\e=eF and is the zero-T dynamic 
Drude conductivity, — e'^K,^{ep)/2'Kh^. Note that 
the only source of nonlinearity of CTpii with respect to the 
ac field power in Eq. (|ll|l is a dependence of Tin on E^j- 

Alternatively, assuming the dominant role of soft in- 
elastic scattering with energy transfers much smaller than 
T, we can write I\n in the Fokker-Planck form: 



2K^ 



£T 



1 + (W±W,)V2 



(6) 



The photoconductivity at Tin ^ r is thus completely 
determined by /, i.e., in this limit the ac field modifies 
the dc current through the heating: 



fph 



2T:h? 



deKodef 



(7) 



The function /(e) changes abruptly around the Fermi 
energy ep on a scale 



where 



re = max{r,A4, (8) 



Aft. ^ {e£^hj2)[K^{ep)/Ko{ep)Y''' (9) 



and lia — VF[Ko{eF)Tia/ (-fY^"^ is the inelastic length 
[vp is the Fermi velocity). Note that r^n in the regime 
of strong heating (A/i ^ T) should be found self- 
consistently with A/j and thus depends on 

Turning to the evaluation of cTph under the assumption 
that Te <C ef , we first notice that a seemingly reasonable 



InJ = de{B[def+T-^f{l-f)]} 



(12) 



where B{e) — (^{5e)'^W{e,5e)) /2 is the diffusion coef- 
ficient in energy space, W is the corresponding rate of 
inelastic processes, and (...) denotes averaging over the 
energy transfer Se. Equation Q becomes then first-order 
in 9e, which gives /(e) described by the Fermi distribu- 
tion with the effective electron temperature 



Teff = T + A 



FP, 



where 



and 



Afp = e'£tTK^{ep)/2mB{eF) 



(13) 



(14) 



2T:h? 6 
e2 ^2 



2T:h? 6 



(2T-|-AFp)AFpXr (15) 



The case of typical energy transfers ~ T may be qualita- 
tively described by either model with _B(eF)'''in ~ Tf.T . 
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The microwave power and temperature dependences of 
CTph can be found from Eqs. ((TT|l . ((T3|) for a variety of scat- 
tering mechanisms. If one assumes that Tin is determined 
by scattering on acoustic phonons via the piezoelectric 
interaction screened by the 2DEG, the characteristic en- 
ergy transfer is Te and r^^ cx Tj* . It follows then from 
Eq. that the heating at Ah ^ T is characterized by 
Te oc By using Eq. ^ we get a^i, - oc S^T-^ 

for Ah ^ T and T independent (Tph — cr^ oc st^'^ other- 
wise. 

Having identified the main contribution to tXph in the 
limit Tin/r 1 [diagrams (a),(b) in Fig. it is instruc- 
tive to compare this contribution with that correspond- 
ing to other diagrams [diagrams (c)-(i)]. While the for- 
mer is related to the heating of electrons by the ac field, 
the latter can be regarded as an effect of radiation on 
the impurity scattering and thus represents a classical 
analog of the quantum effect considered in Refs. IgHTll 
Following the procedure given by Eqs. l(2Jl, © and mak- 
ing use of the explicit matrix form of the field operator. 



nn' 



= ej, + i{v' - v)ey , 



(16) 



one can readily calculate the photoconductivity at any 
desirable order in the fields E^c and £^ {'^x,y are the unit 
vectors along the x, y axes). At the lowest order the pho- 
toconductivity CTph 0{£'^^£'^) is given by the diagrams 
(a)-(i) in Fig. Q (together with their counterparts in a 
lower half-plane, uj — *■ —uj). The result takes a simple 
form in the limit rr} <C r^^ ^ ^ 



Cph 



2ci- 



1 D / eE^vp 

£pUJ 



(5ci -I- 4c2) + (3ci + 2c2) COS 2(/)is 



(17) 



where ci = erd'^eT ^je^jj^, C2 = erd^T '^\e=eF are num- 
bers (typically of order unity) determined by the type of 
disorder, (jyp is the angle between S^c and S^. The first 
term in the square brackets corresponds to the diagrams 
(a), (b) in Fig. ^ and reproduces Eq. (fTT|l in the hmit of 
weak heating. Ah ^ T. The term (5ci-|-4c2) corresponds 
to the diagrams (c), (e), and (f). The polarization- 
dependent part, given by the last term, originates from 
the diagrams (d) and (i) (in which both diagonal links 
have the same direction along the v axis). Finally, the 
diagrams (g) and (h) give a contribution which is smaller, 
compared to the diagrams (c)-(f) and (i), in the param- 
eter l/wcTin ^ 1 and is omitted in Eq. ((T^ . One can 
clearly see from Eq. H18|) that in the limit Tin/r ^ 1 the 
photoconductivity is dominated by the heating of elec- 
trons. 



III. INTERACTION-INDUCED 
PHOTOCONDUCTIVITY 

In the above, we have neglected inelastic electron- 
electron collisions, whose role is not essential for the qua- 
siclassical OPC, but have also ignored the renormaliza- 
tion of the elastic scattering rate by electron-electron in- 
teractions. The latter approximation, which fits in with 
the conventional approach to the photoconductivity, in 
fact misses an important contribution to CTph. Recall that 
the change of the conductivity due to radiation at Tin ^ r 
comes mainly from the heating. It is most illuminating 
to focus on the model of Eqs. H12|) . H15|I . within which 
Cph — erf is simply proportional to T^^ — . Clearly, 
this contribution to CTph is associated with the term in 
the Drudc conductivity that is quadratic in the small pa- 
rameter T/ep. Substituting T^s for T in the Drude term 
yields CTph given by Eq. (|15|l . On the other hand, there 
are T dependent quantum corrections to the conductiv- 
ity, neglected above, in which one should similarly change 
T — > Teff. At low T, the terms in CTph coming from these 
quantum corrections may easily become larger than the 
classical contribution (|15|l . as we now demonstrate. 

For high-mobility samples, we are mostly interested in 
CTph at not too low temperatures Tt /h ^ 1. In this "bal- 
listic" regime, the most important T dependent term in 
the conductivity at zero B, for the limiting case of short- 
range disorder (r ts), is related to screening of the 
disorder by Friedel oscillations, which translates into a T 
and e dependent renormalization of the elastic scattering 
rate. This quantum interaction-induced term is given 
by ACTint — aie^ I'kY^YTts^ Here a is the interaction 
coupling constant, equal to unity for the Coulomb inter- 
action (under the assumption that k~^^ is much smaller 
than the static screening length). Remarkably, ACTint is 
linear in T/ep, in contrast to the classical T dependent 
term which is quadratic in T jtp. Substitution T — > Toff 
in ACTint yields an interaction-induced term in the pho- 
toconductivity, ACTph ~ aieP' / ■nh?){Tcfi — T)ts, where 
Toff — T = App is given by Eq. ((1^ . For finite B, assum- 
ing that T 3> huJc, fi/r, this term in CTph reads 



Act 



ph 



a — (TcfT 



T)ts 



1 



2^2 



c'S 



(1+^^1)2 



(18) 



This result is obtained by inverting the resistivity tensor 
for which the leading (for T 3> fi-Wc, fi/T") interaction- 
induced correction to pxx is B independent, while that 
to p.j.y may be neglected. It follows from the comparison 
of Eqs. H15(I . H18|I that this quantum contribution to the 
photoconductivity is much larger than the classical one 
provided the effective temperature is low. Toff <ti aep, 
which is satisfied for a ~ 1 in the whole range of temper- 
atures in a degenerate Fermi system. Thus sufficiently 
strong interactions have the effect of greatly enhancing 
the photoconductivity. 

For stronger magnetic fields, fkUc ^ T, another 
mechanism of the interaction-induced photoconductiv- 
ity becomes relevant, related to the interplay^^ of 



5 



quasiclassical memory effects and electron-electron in- 
teractions. For the two-component model of disor- 
der, assuming, as above, that rs ^ r^, the T de- 
pendent correction to the conductivity is AtTint ~ 
a(eVn)(TL/Ts)i/2(TTs/fi)-i/2.i7 with numerical fac- 
tors included, this yields a contribution to the photo- 
conductivity 



Act 



ph 



-a 



3C(3/2) /r^y/' T-'^^ - T-^^ 



(19) 



Comparing Eqs. ifT^ and (fH^ . one sees that the lat- 
ter mechanism gives a larger contribution to Cph in the 
whole temperature range T < fiwc- At T ~ fiwc, the 
term (|19|) is still larger than that given by Eq. p8|l by 
a factor (wc^l)^^^ S> 1. With increasing T, however, 
the memory-effects induced correction falls off rapidly, as 
exp[—4:TT'^T/ujc], so that at T ~ hujc\n{ujcT]^) a crossover 
to Eq. ifTFI) occurs. 



IV. MAGNETOOSCILLATIONS DUE TO 
MEMORY EFFECTS 

The photoconductivity obtained in Sees. IIIIlTll exhibits 
the cyclotron resonance but shows no oscillations with 
varying u/ujc- Let us now incorporate the memory ef- 
fects. To this end, we have to step back to write the 
Liouville equation not yet averaged over positions of im- 
purities. More precisely, for the two-component model of 
disorder (specified in SecP), we average only over smooth 
disorder and represent the Liouville operator as 



L = Lo + 6L — /j, 



(20) 



where Lq includes the effect of scattering on smooth dis- 
order: 



Lo = dt+ vVr 



(21) 



and 5L — — J^i -^R; (i") is a sum of collision operators for 
short-range impurities located at points R^. We have to 
keep in Lq the spatial-gradient term (v is the velocity). 

Averaging the solution of Eq. Q over with L given 
by Eg. II20II can be done systematically along the lines of 
Ref. 13 a classical diagram technique is formulated by 
means of the free propagator Lq^ and the disorder corre- 
lation function {SL{r)SL{r')) . We proceed by represent- 
ing the averaged propagator {L^^) — (Lo + -M — ^in)~^ in 
terms of the self-energy operator M. Equations ©-(CI) 
are then reproduced with r in Eq. © given by 



.-1 _-i 



(22) 



where S = J{d(p/2TT) nMn and n = v/|v|. To first order 
in SL, 



-(1) _ 



-ns 



dr J (d0/2^)n/R,(r)n . 



(23) 



By definition ns J dr/R. (r)n = —riTg^, so that we have 
S^^) = Tg^, which yields the Drude result for the total 
scattering rate r^^ = + Tg^. Expanding now M to 
second order in 5L, we obtain the leading correction to 
S that is due to the memory effects: 

= -nsjdrjdr'j^ 
X n/R,(r)i?^(r-r')/R,(r')n , (24) 
or more explicitly 

X cos(j)I-R,i{r,(j),(l)) D^{v - t' ,(j),(j}') 

X /R.(r',^',0') COS0' , (25) 

where the propagator 



(26) 



is taken in the to representation. Most importantly, the 
w dispersion of D leads to oscillations of Yi^'^\u!) with a 
period ujc- 

To find S(2)(w), we first note that, since /Ri(r) as 
a function of r falls off fast beyond a small vicinity 
of Rj, one can put r = r' in the argument of D,^ in 
Eqs. (I11I),(|23. Then, S(2)(tj) is given by 



tistJ J 27r^ 
X cos cj)D^[0;(j),(l>') cos (j)' , (27) 

where Di^{0; (p, (j)') is the Fourier transform in t of the 
probability density to return with a direction of v speci- 
fied by (f)' if one starts at an angle (j). Let us now focus on 
the case uJcTs ^ 1. In this hmit, Dui{0;(j),(l)') is sharply 
peaked dX 4> — 4>' and, introducing the total probability 
of return Pui = J dcp' D^{Q] (j), (/>'), we finally get 



(28) 



A return-induced correction to the effective scatter- 
ing rate, which comes according to Eq. H28|l from RePc^, 
yields, away from the cyclotron resonance, oscillations of 
the absorption rate through a correction to the function 
RUep) [cf. Eq. ©]: 

AK^iep) = -ieF/2nsTl)ReP^/iu; ±u;,)^ . (29) 

The oscillatory part of K^^ep) leads to classical oscilla- 
tions of the linear ac conductivityji^ 



(30) 



and being substituted in Eqs. lO-Q, to those of (Xph. 
To first order in £^ , the classical oscillatory correction to 
CTph reads 



^ph 



Act, 



ReP^ 
nsTs 



(31) 
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(a) 



(b) 



(L )„, 



(l')oo (l')„ D(q,CO) {l\ 



density Pn{x±,x\\) for particles on the Fermi surface to 
be scattered from the starting point on a cyclotron orbit 
by a distance x± across the orbit and a distance X|| along 
it in time 27Tn/uJc is then given by the anisotropic Gaus- 
sian distribution with averages {x\) = /3 — 

(see Appendix). Summing over multiple cyclotron revo- 
lutions, wc thus express P^^ as 



P.. 



E 

ji=i 



(C) 



<- 



( L ),o D(q, 0) D(q, 0) D(q, CO) ( L 



FIG. 2: Diagrams describing the memory effects in the pho- 
toconductivity (jph: oscillatory self-energy (b) and vertex (c) 
corrections to the smooth part (a) of CTph. 



It is worth noting once more that both the smooth cor- 
rection CTph — ctJ* and the oscillatory contribution Aa'^^ 
are proportional to the inelastic time Tin. 

In the above, we have analyzed the oscillatory correc- 
tion to the self-energy in terms of the return probability 
Pui- In fact, there are other contributions to the OPC 
which are not reduced to the self-energy corrections and 
cannot be represented through P^. To illustrate this 
point, it is convenient to switch to a more conventional 
(dual) representation of the diagrams in Fig. ^ now with 
lines corresponding to the propagators and vertices rep- 
resenting the field operators as shown in Fig. |21 
The diagram (a) in Fig. |21 reproduces the graph (a) in 
Fig. n The diagram Fig. |2Ib) represents the oscillatory 
correction to CTph of the self-energy type, Eq. Both 
diagrams (a) and (b) in Fig. |21 contain the inelastic prop- 
agator (L~^)oo = Tin at zero momentum which is much 
larger than all other propagators, {L~^)i,n with at least 
one of the indices v,n ^ 0. By contrast, the diagram 
(c) , which exemplifies an oscillatory vertex correction to 
CTph, is not proportional to Tin, because of large q run- 
ning along the internal propagators D(q,uj) [defined in 
Eq. H26|) ]. The vertex type corrections, which are of the 
same order in all of the diagrams (a)-(i), are thus by 
a factor Tin/rg smaller than the self-energy contribution 

®- . 

The function for lJcTl 3> 1 is most directly evalu- 
ated by using Eq. ipT|l which represents the time evolu- 
tion of (/) as a diffusion process with a white noise spec- 
trum of dfcj). This approach is justified for not too strong 
B, namely for 5 ^ d, where (see Appendix) 



5 = 27ri/2„^ri/(c^,Ti)3/2 



(32) 



is a mean-square fluctuation of the guiding center of a 
cyclotron orbit after one cyclotron revolution (otherwise 
adiabatic drift dynamics is developed). The probability 



Note that once the particle hits a short-range impurity, 
its guiding center is shifted by a distance of the order 
of the cyclotron radius. As a result, the contribution 
of such trajectories to the return probability can be ne- 
glected and only non-colliding orbits should be taken into 
account, which is expressed by the exponential factor 
exp(— t/rs). Equation gives oscillations of P^ as 
Lu/ujc is varied: 



P, 



1 



1 



^ n — 1 



: exp 



(ito + F) 



(34) 



whose damping with decreasing B is characterized by 

r = ^{ — \ +^ . (35) 



_3_ 
2^ 



OJr 



In the limit of weak damping, ttF <C lOc, we perform 
the summation in Eq. (|34|l by means of Poisson's formula 
to represent Re P^^ as a series of sharp peaks centered at 
u! — NlUc- a peak at ~ Nuic '° ^^""^ f"'-"^ 



RePc^ 



2\/3Trv pU! 
1 + 



is of the forma* 



1 



.2^1/2 



2(1 



F 

1/2 



(36) 



T{0) = 1 . (37) 



Note that the amplitude of the peaks in Eq. 136|) falls off 
with decreasing lUc or increasing u; as a power law, namely 
as i^c/^- The power-law suppression of the oscillations 
crosses over into the exponential damping only for very 
large Tr(uj/ujc)'^ ^ ^c^l, when one can neglect all terms 
m Eq. (PU but the first one, which gives 



ReP^ = 



[u^TLf/^ 2TTUJ 

5 cos exp 

2'KVpTL LOc 



27rF\ 

) 



(38) 



It is worth noting that, because of the condition uJcTs 2> 
1, the term Tg^ in Eq. (|35|) may be neglected in the damp- 
ing factor of Eq. I|38|l . so that the exponential damping is 
determined by the momentum relaxation time for scat- 
tering off the long-range disorder. 



V. OSCILLATORY AC CONDUCTIVITY: 
QUASICLASSICAL VS QUANTUM 

Now we compare the classical oscillatory ac conductiv- 
ity ffif^ given by Eqs. (|Sn|l . . . with the quantum 
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%^ 0,8 




FIG. 3: Quasiclassical [o-^f', Eq. 1401 1 and quantum [a-i?\ 
Eq. (I4lll 1 oscillatory ac conductivity (normalized to the Drude 
conductivity a^) vs luc/lo for uj/2tv — 100 GHz, r = 0.6 ns, 
r/rq = 50, ts/tl = 0.1, a/S = 0.25 at tJc/t^ = 1/2. 



contribution tr, 



(9) 



(TiL"' for weak damping at oj 



calculated in Ref. 7. Let us represent 

Nur as 



D 



1 



1/2 



(39) 



where (5 is given by Eq. (|32(l and we have introduced 
a characteristic size of the short-range impurities a = 
{nsVFTs)~^- It follows that, apart from the harmon- 
ics number N, the amplitude of the oscillations is given 
by the product of a small factor a/ 6 and a large factor 
(uJcTl)^/'^. In the exponential damping regime, a!^"^ 
re-written as 



IS 



= 1 - 



2nu} 
cos exp 



37r 

l^cTL 



, (40) 



so that the pre-exponential factor is simply given by a/ 6. 
An important point to notice is that the damping in 
Eq. (|40|l is characterized solely by the long transport 
time for scattering off the smooth disorder. On the other 
hand, the envelope of the quantum oscillations of the ac 
conductivity is determined by the single-particle time r^: 



1 



2t:uj 
2 cos exp 



2tt 



(41) 



(this equation is valid for 2ttT ^ fi/Tq, for smaller T see 
Ref.0). 

Note the difference in the sign of the oscillatory terms: 
there is a tt shift of the quantum and classical oscillations 
with respect to each other. Another difference is that the 
damping of the classical oscillations is lu dependent, in 
contrast to the quantum case. One sees that, despite the 
small factor a/S, the classical oscillations may be stronger 
than the quantum ones since in high-mobility structures 
Tq <C tl and the quantum oscillations are damped much 
more strongly. The behavior of the two contributions to 
the oscillatory ac conductivity is illustrated in Fig. |51 



VI. MECHANISMS OF THE OSCILLATORY 
PHOTOCONDUCTIVITY: QUASICLASSICAL VS 
QUANTUM 

Having found the classical contribution Affp'^^^ to the 
OPC [Eqs. (EU, (ESI), ESI], let us compare it with the 
quantum oscillatory contribution Adp^' related to the 
oscillations of the DOS. Using Eqs. ED , (EH , EHJ and 

omitting numerical factors, we write down the essential 

(c) 

factors in Acpj^ for ttT > ojc for the case of noninteracting 
electrons: 



X cos 



The sign of ACTpj^ in Eq. (gSJ depends on that of Kq [see 
Eq. ini)]- 

As shown in Sec. lIIII unless the electron-electron inter- 
action is extremely weak, the largest contribution to the 
smooth part of CTph comes from the interaction correc- 
tion to the conductivity. For T > f>iOc^'n{uJcTL), the main 
interaction-induced term in the smooth part of Cph is 
given by Eq. (fl^ and, according to Eq. EU, this yields in 

(c) 

turn the main term in the oscillating part AcTpj; . Putting 
a ^ I (long-range Coulomb interaction), we have: 







2 








2ttuj 










exp 




(- 











(42) 



Ad. 



ph 



D 'in 



eFT 



2HL0 

cos exp 



37r 



(43) 



which is larger than the noninteracting part [Eq. H42|) ] by 
a factor ep/T. 

For T <C hujc, the main contribution to (Tph is re- 
lated to the interaction correction given by Eq. H19|l. 

which yields the oscillatory part AUp'jj' similar (in terms 
of the phase of the oscillations and their damping fac- 
tor) to that in Eq. (|43|l but multiplied by a large factor 
{ujcTlY / (Ttl / hY^"^ . In the intermediate range of tem- 
perature, fuvc <C T ^ /iWc In (wcTl), there is an expo- 
nentially fast crossover between the two regimes. The 
regime most relevant to the experimentaiiSiMi^ is that 
of high temperature, T > fiWc- It is worth noting that 
Eqs. 142|) . H43|I remain valid in the regime of strong heat- 
ing as well, provided the effective electron temperature 
Te [Eq. ®] is substituted for T. 



For convenience, we also reproduce here Act 



(9) 
ph 



the case of overlapping Landau levels; specifically, for 
|w ± Wcj ^ LOc in the regime linear with respect to the 
microwave power (see Eq. (8) of Ref. 8; here we omit 
numerical factors): 



Aa' 



(9) 
ph 



D 'cc 



cE^vf 
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X sm exp 



27r 



(44) 



The electron-electron scattering time^ t^c oc (up to 
a logarithmic factor) depends on the effective electron 
temperature [Eq. ©]. Although both contributions, 
Eqs. (|43|l ■ (|44(l have the same period in uj/uJc, crucial dis- 
tinctions are clear. 

Firstly, their phases are shifted by n/2. Secondly, 
despite both contributions being proportional to a cer- 
tain inelastic relaxation time, they are different in that 
the amplitude of Aa^^ is limited by t^q (which at low 
T is much shorter than the electron-phonon scattering 
time) , whereas the classical term is not sensitive to the in- 
elastic electron-electron scattering in any essential wajii^ 
and is proportional to the energy relaxation time [tih in 
Eqs. ifTTiji . . limited by couphng to the external bath 
(phonons). It follows that in the limit of small T the ratio 
of the amplitudes of the OPC, classical-to-quantum, con- 
tains a large T dependent factor Tin/Too, which may be 
easily as large as 10^. The sensitivity of Aa'^^^ to electron- 
electron collisions stems from the fact that the quantum 
contribution is due to a radiation-induced change of the 
distribution function /(e) that oscillates with both e and 

LV. By contrast, the classical contribution l^a'^^ is asso- 
ciated with an oscillatory term in the characteristic elec- 
tron temperature, i.e., with a smooth part of /(e) which 
oscillates with w onlyiSi 

Thirdly, the dependences of the envelope of the OPC 
on uj, bJc, and the degree of disorder are quite different. 
The most important point is that although there is a 
small factor c>c ep^ in Eq. (|43() . in addition to another 
small factor a/5, the damping of the classical term is 

much weaker than that of AtTp''' oc exp(— 27r/ct;c'''q). In- 
deed, the exponential damping of Acpj^^ is governed by 
tl [Eq. (|43() ]. which is far larger than Tq in high- mobility 
samples. It is only that in the limit of very low B that 
the factor in the exponent of Eq. suppresses 
the classical OPC more effectively than the linear in lo'^^ 
Dingle factor in the quantum case. 

It is important to stress that the amplitude of the 
classical OPC in units of the dark conductivity is not 
large under the conditions of the experiments on the 
zero-resistance states. Indeed, the pre-exponential fac- 
tor of Eq. (gnil may be written as eFT){a/ 5) for 
the regime linear with respect to the microwave power. 
Now, the crossover to the regime of strong heating oc- 
curs when the classical OPC is still small, namely the 
ratio Acr^^j^Vc^? is of order {T/eF){a/5) < 10'"^. This 
should be contrasted with the quantum OPC which may 
become large (and thus lead to the zero-resistance states) 
when the heating may be still negligible. For the regime 
of strong heating, when the effective electron tempera- 
ture Te > T, the amplitude of the classical OPC shows 
a sub-linear growth with increasing microwave power 
and may be estimated as a^{Te/eF)ia/S). In particu- 
lar, for the piezoelectric mechanism of the energy relax- 



ation due to electron-phonon coupling, the classical OPC 
grows as £^^^ [see the discussion below Eq. ((T^ ]. We 
conclude that, because of the slow growth with increas- 
ing microwave power, the characteristic ratio Aa^^^/a^ 
can hardly exceed the level of a few percent in the cur- 
rent experiments. That is to say the zero-resistance 
states are related to the quantum OPC. The most favor- 
able conditions for the observation of the classical OPC 
should be realized if the quantum contribution to the 
OPC is strongly damped, which means sufficiently large 
2TT/LUTq > 7. For a high-mobility sample with Tq ^ 3ps, 
this would require uj/2tt < 50 GHz. 

The above analysis shows that the classical OPC can- 
not possibly explain the experimentally reported strong 
deviations of the ratio Tph/rq from 1/2, the value pre- 
dicted by the theory of the quantum OPC, as discussed 
in Introduction. We thus argue that the experiments 
on the damping of Shubnikov-de Haas oscillations might 
strongly overestimate the single-particle scattering rate 
Tq^. One of the reasons could be the presence of 
macroscopic inhomogeneities leading to an inhomoge- 
neous broadening of Landau levels, which might be by 
far larger than the homogeneous broadening given by 
Tq^ and measured in the photoconductivity experiments 
(such a possibility was mentioned in Ref. Jj. We sug- 
gest that measuring the damping of the OPC provides a 
reliable means of extracting from the magnetooscil- 
lations, free from the effect of the additional inhomoge- 
neous damping characteristic to the Shubnikov-de Haas 
measurements. The method based on the OPC is par- 
ticularly useful in high-mobility samples, where is 
small and in the conventional Shubnikov-de Haas mea- 
surements one has to go to fairly low temperatures to 
separate the impurity-induced damping from that related 
to the thermal smearing of the Fermi surface. 



VII. CONCLUSIONS 

In summary, we have analyzed the quasiclassical mech- 
anism of magnetooscillations in the ac- and photoconduc- 
tivity, related to non-Markovian dynamics of disorder- 
induced scattering of electrons in high-mobility struc- 
tures. We have calculated the leading contribution as- 
sociated with a radiation-induced change of the electron 
distribution function, which is proportional to the inelas- 
tic (electron-phonon) relaxation time. We have found 
that the quasiclassical oscillations in the photoconductiv- 
ity are weak under the conditions of current experiments. 
Therefore, the zero-resistance states and the strong os- 
cillations that have been observed in the experiments are 
likely due to the quantum mechanism of Refs. '^'3. We 
argue that the damping of the oscillatory photoconduc- 
tivity provides a reliable method of measuring the homo- 
geneous broadening of Landau levels (single-particle scat- 
tering rate ) in high- mobility structures (which also 
resolves the dilemma posed in Introduction: the anal- 
ysis of the damping of Shubnikov-de Haas oscillations 
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apparently gives overestimated values of due to an 
inhomogcneous broadening) . 

On the other hand, we have identified a range of 
parameters within which the quasiclassical mechanism 
yields oscillations of the photoconductivity that may 
dominate at small B over those based on the Landau 
quantization. In addition to the different low-_B damp- 
ing factor, the quasiclassical oscillations are shifted in 
phase by 7r/2 with respect to the quantum oscillations, 
see Eqs. We have also shown that the qua- 

siclassical magnetooscillations in the ac conductivity are 
much stronger than in the photoconductivity and may 
easily compete with the quantum oscillations Z 
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thank I. L. Aleiner for attracting our attention to the im- 
portance of quantum interaction corrections. This work 
was supported by the SPP "Quanten-Hall-Systeme" of 
DFG and by RFBR. 



APPENDIX: RETURN PROBABILITY IN A 
MAGNETIC FIELD 

The return probability [Eq. (|28|) ] can be directly 
evaluated by using the quasiclassical propagator D = 
(Lq + '''s ^)~^ [^'l- 03]- Appendix, we present a 

different, more illustrative way to derive P^^. We recall 
that the Liouville operator Lg [Eq. (PU] represents the 
time evolution of the direction of the electron velocity 
n — v/|v| — {— sini/), cos0) as a combination of the cy- 
clotron motion and the angle diffusion due to scattering 
off smooth disorder. The random part x of the angle 
is characterized by a white noise spectrum of dtX'- 



(l){t) = 00 + io^t + x{t) , 
{dtx{t)dt'x{t'))^-5{t-t'). 

TL 



(A.l) 



In what follows we calculate the mean-square fluctuation 
of the guiding center of cyclotron motion S, Eq. (|32() . and 
mean-square shifts of an electron along and across the 
cyclotron orbit after n cyclotron periods at t = nTc — 
27r n/wc- For definiteness, let the guiding center be ini- 
tially placed at the origin, R(< = 0) = (0,0), and the 
electron coordinate and velocity be r (i = 0) = (i?c, 0), 
v{t = 0) = (0, vp) {Rc = vp/uJc is the cyclotron radius). 



Using Eq. (jA.ip we get the mean-square shifts of the po- 
sition of guiding center R and the fluctuating angle x in 
time t = nTr'. 



{Rl) = 




TL Jo 

{Rl) = {Rl) , 



dtRc cos (t>{t)dtx{t) 



^"'^diCOs20(i) 

TL 




Atdtxit) 



2nTc 

TL 



(A.2) 



The diffusion approximation is valid as long as the root- 
mean-square shift of the guiding center after one cy- 
clotron revolution (n = 1) exceeds the characteristic 
length scale of the random potential. 



5^[{Rl) + {Rl) 



^^^"'-RA^\'\>d. (A.3) 
^cTl J 



In the same manner, we calculate the mean square of 
electron shifts along and across the cyclotron orbit, = 

y{t = nTc) = Jq'^" dtvp cos(/)(t) and xj^ = x{t = nTc) — 

Rc~— Jo'^" dtvp sin(/)(t), respectively. 




dtvp cos{LUct) / dt' dt'xit') 
Jo 



{xp = {{Ry + i?, xY) = {Ri) + Ri {x') ■ (A.4) 

It follows that fluctuations along the cyclotron orbit 
are enhanced with respect to those across the orbit, 
(a;|) = 3(a;^) = 3 n 5^/2, and we arrive at the anisotropic 
electron distribution after n cyclotron revolutions. 



Pn{x±,Xn) 



1 



1 



exp 



: exp 



nS 
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which enters Eq. H33|l for the return probability 



, (A.5) 
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More precisely, the effect of inelastic electron-electron scat- 
tering on the OPC is twofold: firstly, it suppresses oscilla- 
tions of f{e) induced by those of the DOS; secondly, it af- 
fects the exact shape of the smooth part of /(e) by thermal- 
izing electrons among themselves. The former effect leads 
to the damping of the quantum term Ao-p*^' , as discussed 
in Ref. 0. The latter can only yield a numerical coefficient 
of order unity in the classical term Aa^^^^, depending on 
what other mechanisms of inelastic scattering are. More- 
over, this numerical factor is one, so that electron-electron 
inelastic scattering does not manifest itself in Aa'^^^ at all, 
if the inelastic coupling to the thermal bath in the pres- 
ence of a driving force yields /(e) whose shape is given by 
the Fermi distribution (but with an effective temperature 
different from that of the bath). This is the case, e.g., for 
the Fokker-Planck mechanism of Eqs. (1121 . 115^ . 
Similar oscillations of the ac conductivity were studied for 
the case of a random antidot array, where the memory ef- 
fects are still stronger, in D.G. Polyakov, F. Evers, and 
LV. Gornyi, Phys. Rev. B 65, 125326 (2002). The antidot- 
array model also describes correctly rare strong short- 
ranged scatterers in the absence of background smooth 
disorder. 

Equation 1361 1 has the same structure as an interaction- 
induced correction to the tunneling DOS, studied in 
A.M. Rudin, LL. Aleiner, and L.I. Glazman, Phys. Rev. 
Lett. 78, 709 (1997), which is also proportional to the re- 
turn probability P„. 

It is worth noting that there exists also a subleading quan- 
tum contribution to the OPC that is due to the oscillatory 
heating, i.e., due to a smooth part of /(e) which oscillates 
with uj (because of the oscillatory absorption rate coming 
in turn from the Landau quantization). It has the same 
damping factor as in Eq. 1441 , but the phase of the oscilla- 
tions as in Eq. 1431 . This quantum term is, however, much 
smaller than given by Eq. 1441 : specifically, by a factor 

{h^UJUJc/e%)Tin/Tcc <. 1. 



